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.^^ ' Abstract 

^ ■ We consider the compressible (barotropic) Navier-Stokes system on time-dependent domains, supplemented 

-4-^ ' with slip boundary conditions. Our approach is based on penalization of the boundary behaviour, viscosity, and 

^ - the pressure in the weak formulation. Global-in-time weak solutions are obtained. 
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^^ 1 Introduction 

, ' Problems involving the motion of solid objects in fluids occur frequently in various applications of continuum fluid 

• . dynamics, where the boundary conditions on the interfaces play a crucial role. Besides the commonly used no-slip 

SIJ ■ condition, where the velocity of the fluid coincides with that of the adjacent solid body, various slip-like conditions 

*vj . have been proposed to handle the situations in which the no-slip scenario fails to produce a correct description of the 

fluid boundary behavior, see Bulicek, Malek and Rajagopal |Tj, Priezjev a Troian [TT] and the references therein. For 

viscous fluids, Navier proposed the boundary conditions in the form 

>< : [Htan + ^l^- V]ta„ lr. = 0, ^ > 0, (1.1) 

. . .' where § is the viscous stress tensor, k represents a "friction" coefficient, u and V denote the fluid and solid body 
velocities, respectively, and Pi- is the position of the interface at a time r, with the outer normal vector n. If k = 0, 
we obtain the complete slip while the asymptotic limit k — > oo gives rise to the standard no-slip boundary conditions. 
Besides their applications in "thin" domains occurring in nanotechnologies (see Qian, Wang and Sheng [12), the 
slip boundary conditions are particularly relevant for dense viscous gases (see Coron [2]), described by means of the 
standard Navier-Stokes system: 

dtg + diY,{gu)^0, (1.2) 
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dtigu) + diva;{gu u) + "^xPig) = div^§(V:ru) + gf, (1.3) 

where g is the density, p = p{g) the (barotropic) pressure, f a given external force, and § is determined by the standard 
Newton rheological law 

S(Vxu) = fi ( VxU + V^u - -divajul J + 77div2;Ul, /i > 0, 77 > 0. (1.4) 

The boundary of the domain fit occupied by the fluid is described by means of a given velocity field V(t, x), where 
t > and X E R^. More specifically, assuming V is regular, we solve the associated system of differential equations 

^X(i,x)=v(t,X(i,x)), t>0, X(0,x)=x, (1.5) 

and set 

fir ~ "X. (r, flo) , where flo C i?^ is a given domain, F^ = dflr, and Qr = {{t, x) | i G (0, t), x <E fir}- 
In addition to (jl.ip . we assume that the boundary Tr is impermeable, meaning 

(u- V) -nlr, =0 for any T> 0. (1.6) 

Finally, the problem (|l.ll - II. 6p is supplemented by the initial conditions 

£^(0, •) = go, {gu){0, •) = {gu)o in Hq. (1.7) 

Our main goal is to show existence of global-in-time weak solutions to problem (jl.ll - I1.7P for any finite energy 
initial data. The existence theory for the barotropic Navier-Stokes system on fixed spatial domains in the framework 
of weak solutions was developed in the seminal work by Lions jT^] , and later extended in [2] to a class of physically 
relevant pressure-density state equations. The investigation of incompressible fluids in time dependent domains started 
with a seminal paper of Ladyzhenskaya [11] , see also [13| , [M] , [15] for more recent results in this direction. 

Compressible fluid flows in time dependent domains, supplemented with the no-slip boundary conditions, were 
examined in [7] by means of Brinkman's penalization method. However, applying a penalization method to the slip 
boundary conditions is more delicate. Unlike for no-slip, where the fluid velocity coincides with the field V outside 
fl-r, it is only its normal component u • n that can be controlled in the case of slip. In particular, given the rather 
poor a priori bounds available in the class of weak solutions, we lose control over the boundary behavior of the normal 
stress §n involved in Navier's condition (jl.ip . 

A rather obvious penalty approach to slip conditions for stationary incompressible fluids was proposed by Stokes 
and Carey |20| . In the present setting, the variational (weak) formulation of the momentum equation is supplemented 
by a singular forcing term 

- / / (u - V) • n V5 • n dS:i, di, £ > smaU, (1.8) 

£ Jo JTt 

penalizing the normal component of the velocity on the boundary of the fluid domain. In the time-dependent geome- 
tries, the penalization can be applied in the interior of a flxed reference domain, however, the resulting limit system 
consists of two fluids separated by impermeable boundary and coupled through the tangential components of normal 
stresses. In such a way, an extra term is produced acting on the fluid by its "complementary" part outside fir- In 
order to eliminate these extra stresses, we use the following three level penalization scheme: 



1. In addition to (|1.8p . we introduce a variable shear viscosity coefficient /.t = /i(j, where ji^^ remains strictly positive 
in the fluid domain Qt but vanishes in the solid domain Q9p as w — > 0. 



2. Similarly to the existence theory developed in [5], we introduce the artificial pressure 

P&{q)^p{q) + 5q'\ /3>2, J>0, 
in the momentum equation (|1.3p . 

3. Keeping e,5, w > fixed, we solve the modified problem in a (bounded) reference domain B <Z R^ chosen in 
such a way that 

i^T C B for any t > 0. 

To this end, we adapt the existence theory for the compressible Navier-Stokes system with variable viscosity 
coefficients developed in [6]. 

4. We take the initial density go vanishing outside Hq, and letting e — > for fixed S, w > we obtain a "two-fluid" 
system, where the density vanishes in the solid part ((0, T) x B)\ Qt of the reference domain. 

5. Letting the viscosity vanish in the solid part, we perform the limit w — > 0, where the extra stresses disappear in 
the limit system. The desired conclusion results from the final limit process (5 — > 0. 

The paper is organized as follows. In Section [2J we introduce all necessary preliminary material including a weak 
formulation of the problem and state the main result. Section [3] is devoted to the penalized problem and to uniform 
bounds and existence of solutions at the starting level of approximations. In Section 31 the singular limits for e -> 0, 
w — > 0, and (5 — > are preformed successively. Section [5] discusses possible extensions and applications of the method. 

2 Preliminaries, weak formulation, main result 

In the weak formulation, it is convenient that the equation of continuity (jl.2p holds in the whole physical space R^ 
provided the density g was extended to be zero outside the fluid domain, specifically 



/ gifi{T,-)dx- / go(p{0,-) dx ^ / / {gdtip + gu ■ W xip) dx dt 

JU-r J^o Jo Jvit 



(2.1) 



for any r S [0,T] and any test function ip S C^([0,T] x R^). Moreover, equation (|1.2p is also satisfied in the sense of 
renormalizcd solutions introduced by DiPerna and Lions |3]: 

b{gMT,-)dx- [ b{goMOr)dx= [ [ {b{g)dt^ + b{g)u ■ V,^ + {b{g) - b'{g)g) div.u^) dx dt (2.2) 

for any r G [0,T], any tp e C'^{[Q,T] x R^), and any b e C^[0,oo), b{0) = 0, b'{r) = for large r. Of course, we 
suppose that g> a. a. in (0,T) x R^. 

Similarly, the momentum equation (|1.3p is replaced by a family of integral identities 

/ gu-^{T,-)dx- f {gu)o-^{0,-)dx (2.3) 

= / {gu ■ dt<fi + g[u (X) u] : V^f + p(£')div^(p — S(Va;u) : VxV + £'f ' v) da; dt 

Jo Jilt 



for any r € [0,T] and any test function ip G C;?°([0,r] x R^;R^) satisfying 

(f ■ n|r, = for any r e [0, T]. (2.4) 

Finally, the impermeability condition (jl.6p is satisfied in the sense of traces, specifically, 

ueL^iO,T;W^'^{R^;R^)) and (u - V) • n(r, •)|r, = for a.a. r G [0,r]. (2.5) 

At this stage, we are ready to state the main result of the present paper: 



Theorem 2.1 Let Jlo CI R^ he a hounded domain of class C^^" , and let V e C^{[0,T];C^{R^;R^)) he given. 
Assume that the pressure p € C[0, oo) n C^{0, oo) satisfies 

p{0) — 0, p{g) > for any g > 0, lim — — r = poo > /or a certain 7 > 3/2. 

Let the initial data satisfy 

Qo e L'^{R^), Qo > 0, go ^ 0, £'oU3\Oq = 0, {qu)q = a.a. on the set {go = 0}, / — |(eu)o|^ dx < 00. 

Jua Qo 

Then the prohlem \L1\ - rOj) admits a weak solution on any time interval (0,T) in the sense specified through i2.1\ 



The rest of the paper is devoted to the proof of Theorem 12. II 

3 Penalization 

For the sake of simplicity, we restrict ourselves to the case k — 0, r] = 0, and f = 0. As we shall see, the main ideas of 
the proof presented below require only straightforward modifications to accommodate the general case. 

3.1 Penalized problem - weak formulation 

Choosing R > Q such that 

V|[o,T]x{|x|>fl}-0, noC{|x|<i?} (3.1) 

we take the reference domain _B = {|x| < 2R}. 

Next, the shear viscosity coefficient fi — ^^{t,x.) is taken such that 

fi^ e C^i[0,T] X i?3), < /^ < fi^{t,x) <^lm [0,T] x B, h^{t, ■)\n^ - fi for any r e [0,r]. (3.2) 

Finally, we define modified initial data so that 

Qo = go,S: go.s >0, go,s^O, go.slB^^XUo =0, / (g2^g + dg^Adx<c, (3.3) 



{gu)Q = {gu)o.s, igu)o^5 = a.a. on the set {go,s = 0}, / \{gu)o,s\'^ dx < c. (3.4) 

Jn„ Qo.s 

The weak formulation of the penalized problem reads as follows: 

/ gip{T,-)dx~ / goip{0,-) dx = / / {gdtip + gu-Vx^) dx dt (3.5) 

Jb Jb Jo J b 

for any r e [0,T] and any test function ^ £ C^{[Q,T] x R^); 

I g\i-ip{T,-) dx- I (giu)o • (^(O,-) dx (3.6) 

JB JB 

= i gu- dtif + g[vi ® u] : V^v' + p{.Q)dYVxV + '^e'^diva;^' — ^l^^ I V^u + V^u — -divj;Ul j : V^t/J j dx dt 

+ - I [ ((V-u)-nip-n) dS^ dt 

for any r G [0, T] and any test function (f e Cc'"([0, T] x B; R^), where u e L^{0, T; Wj '^(B; i?^)), meaning u satisfies 
the no-slip boundary condition 

u\gB = in the sense of traces. (3.7) 

Here, e, S, and lo are positive parameters. The choice of the no-slip boundary condition p.7p is not essential. 

The existence of global-in-time solutions to the penalized problem can be shown by means of the method developed 
in [6] to handle the nonconstant viscosity coefficients. Indeed, for e > fixed, the extra penalty term in (13. 6p can be 
treated as a "compact" perturbation. In addition, solutions can be constructed satisfying the energy inequality 



7;Qm 



PiQ) 



^ gAir,-)dx+l 



(3-1 



JB 



Vxu + V^u - -div^ul 



da; dt 



(3.8) 



~e£ i [("-V).n]u.ndS, di<^(^^|(eu)o,5|2 + P(go,,) + ^e^,j dx. 



where 



Pig) = gJ' 



' p{z) 



dz. 



Note that the quantity on the right-hand side of (|3.8p representing the total energy of the system is finite because of 
dSD, ([331). 

In addition, since /3 > 2, the density is square integrable and we may use the regularization technique of DiPerna 
and Lions [3] to deduce the renormalized version of <\3.5\i . namely 

b{g)ip{T, ■)dx- [ bigoMO, ■) dx = [ [ {b{g)dtV + Kq)^ ' ^x^ + {b{g) - b'{g)g) div.u^) dx dt (3.9) 

B JB JO JB 



for any ip and b as in (|2.2 



3.2 Modified energy inequality and uniform bounds 

Since the vector field V vanishes on the boundary OB, it may be used as a test function in 
resulting expression with the energy inequality (|3.8p . we obtain 



Combining the 



+ ^^')^^'-)'^nioL^^ 



V^u + V^u - -div::cul 



da; dt 



(3.10) 



+ I ((gu • V)(t, •) - {gu)oj ■ V(0, •)) d.T 

+ / / (a^cu fV:rU + V^u--div^uIj -.V^V- gu-dtV- gu(Su:V^Y~p{g)dw^V- -g^TY/div^Vj da; dt 
Since the vector field V is regular and since 

p{g) < c(l + P{g)) for aU g>0, 
relation (j3.10p gives rise to the following bounds independent of the parameters e, i5, and to: 



ess sup ||Veu(^r)llL2(i3;i?3) < c, 

te(o,T) 
ess sup / P(£»)(t, •) da; < c yielding ess sup ||e(i, OIlL-'fs) — "^j 



1/9 



ess sup S\\g{t,-)\\lf,,g. < c, 

tG(0,T) 



JB 



Va;U + Vlu - -div^ul 



da; dt < c, 



and 



/ / |(u-V)-n|^ dS:, di<£c. 
Finally, we note that the total mass is conserved, meaning 

/ g{T, •) da; = / go.s dx = go.s da; < c for any r G [0,r]. 
Jb Jb Juf, 



(3.11) 

(3.12) 
(3.13) 

(3.14) 
(3.15) 

(3.16) 



Thus, relations p. lip . p.l4p . p.l6p . combined with the generalized version of Korn's inequality (see [8j Theorem 
10.17]), imply that 



|u(t,-)ll 



<-^(i3;fl3) 



< c{uj). 



(3.17) 



3.3 Pressure estimates 

Since the surfaces Ft are determined a priori, we can use the technique based on the so-called Bogovskii operator to 
deduce the uniform bounds 

/ / (p{q)q'' + Sg/^+A dx di < c{JC) for a certain i^ > (3.18) 

for any compact JC C [0,T] x B such that 

ICn{Ure[o,T]{{r}xrr)) =0, 

see [TU] for details. 

Note that due to the fact that the boundaries r,- change with time, uniform estimates like p.lSp on the whole 
space time cylinder (0,T) x B seem to be a delicate matter. On the other hand, the mere equi-integrahility of the 
pressure could be shown by the method based on special test functions used in [7]. 

4 Singular limits 

In this section, we perform successively the singular limits e — >■ 0, w — >■ 0, and 5 — >■ 0. 

4.1 Penalization limit 

Keeping the parameters 5, ui fixed, our goal is to let e — >• in p.Sp . (j3.6p . Let {g^,u^} be the corresponding weak 
solution of the perturbed problem constructed in the previous section. To begin, the estimates p.l2p . p.l7p . combined 
with the equation of continuity p.5p , imply that 



&^ginCweak([0,T];i^(B)), 

and 

Ue -> u weakly in L^{0,T;Wq'^{B, R^)) 

at least for suitable subsequences, where, as a direct consequence of ()3.15p . 

(u-V)-n(T,-)|r. =Ofora.a. re [0,r]. (4.1) 

Consequently, in accordance with (|3.1ip . (|3.12p and the compact embedding U'{B) ^->^-> W^^'^{B), we obtain 

Q^xi^ -^ £.u weakly-(*) in L°"{0,T; L^^/'-^+'^\B; R^)), (4.2) 

and, thanks to the embedding W^-^{B) ^ L^{B), 

Q,u, (g) Ue ^ eu(g)u weakly in L^{0, T; l67/(47+3)(5. ^3^^^ 

where we have used the bar to denote a weak limit of a composed function. 
Finally, it follows from the momentum equation p.6p that 

PeUe -> gu 



inC7„cak([ri,T2];L27/(7+i)(0;i?3)) 



for any space-time cylinder 

(ri,T2) X O c [0,r] X B, [Ti,T2] X on U,e[o,T] ({r} x T,) 
Seeing that L27/(7+i)(5) ^^ W^-i'2(B) we conclude, exactly as in (g^l), that 



pu (g) u = gu (g) u a. a. in (0, T) x B. 
Passing to the limit in (13. 5p we obtain 



/ Qip{T,-)(lx- I go,5'p(0, •) da; = / / {gOti^ + gxi ■ \J x'^) dx At 
JB Jb Jo Jb 



(4.3) 



for any r G [0,T] and any test function ip e C^{[0,T] x i?^). 

The limit in the momentum equation (j3.6p is more delicate. Since we have at hand only the local estimates (I3.18P 
on the pressure, we have to restrict ourselves to the class of test functions 

¥>eCi([0,T];<'°°(S;i?3))^ supp[div,<^(T, •)] n T, = 0, <^ • n|r, = for aU r e [0,T]. (4.4) 

Passing to the limit in (13. 6p , we obtain 

/ gu ■ ip{T, ■) dx - / (£»u)o,5 • (^(O,-) dx (4.5) 

Jb J b 

= i gu- dt^ + g[vL (g) u] : Vx^ + p{g)diYxip + 5gf^<:\\Vx^ — ^^ i Va;U + V^u — -divj;Ul J : VxV ) dx dt 



for any test function ip as (|4.4p . Note that the requirement of smoothness of ip postulated in p.6p can be easily relaxed 
by means of a density argument. 

Finally, we show pointwise (a. a.) convergence of the sequence {ge}e>o- To this end, we adopt the method developed 
in [6] to accommodate the variable viscosity coefficient ^^. The crucial observation is the effective viscous pressure 
identity that can be established exactly as in 6j: 



P5ig)Tk{g) - ps{g) Tk{g) = -l^u \Tk{g)d\VxVi - Tk{g)d\v x'aj (4.6) 

where we have denoted 

Ps{q) = p{q) + (5/, Tk{g) = min{gi, k}. 

Similarly to the pressure estimates (|3.18p . identity (|4.6p holds only on compact sets /C C [0,T] x i? satisfying 



/cn(u,e[o,T](Mxr, 



We recall that this step leans on the satisfaction of the renormalized equation (|3.9p for both g^ and the limit g that 
can be shown by the regularization procedure of DiPerna and Lions [3] . 
Following [6j, we introduce the oscillations defect measure 



osCg[£<£ -^ g]{IC) = sup ( limsup / \Tk{gir) - Tkig)]"^ dx dt 

k>0 \ e^O Jk. 



and use (|4.6p to conclude that 

osc^ +i[g^ ^ q]{IC) < c{u!) < oo, (4.7) 

where the constant c is independent of /C. Thus 

osc^+i[e, -> e]([0,r] xB)< ciuj), (4.8) 

which imphes, by virtue of the procedure developed in [5], the desired conclusion 

g, -> g&.a. in (0,r) x B. (4.9) 

In accordance with (|4.9p . the momentum equation (|4.5|) reads 



/ gu • <f{r, ■)dx- f igu)o,s ■ ¥^(0, •) dx (4.10) 

jb Jb 

= / ( ^" ' ^tV' + ^[" ® "] ■ ^^'^ + p{g)dWx<fi + Sg'^div^ip - ^„ f V^u + V^u - -divj;Ul j : Vx<p j dx dt 

for any test function tp as in ()4.4|) . In addition, as already observed, the limit solution {g, u} satisfies also the 
renormalized equation p.9p . 



4.1.1 Fundamental lemma 

Our next goal is to use the specific choice of the initial data go^s to get rid of the density-dependent terms in (|4.10p 
supported by the "solid" part ((0,T) x B)\ Qt- To this end, we show the following result, rather obvious for regular 
solutions but a bit more delicate in the weak framework, that may be of independent interest. 



Lemma 4.1 Let g e L°°{Q,T;L'^{B)), g>Q,u€ 
continuity, specifically, 


L2(0,r;<'2(i?;i?3)) be 


a weak 


■ solution 


of 


the 


equation of 






IM 


,-)(piT,-) -£.0V(0, •))( 


ix ^ / V^^*'^ + ^" ' 


V,(/p) ( 


ix dt 






(4.11) 


for any r € 


[o,r] 


and any test 


function ip G C^([0,T 


'] X i?3). 












In addition. 


assume that 


(u-V)(r,.)-n|r 


V =0 for a.a. re (0,r), 










(4.12) 


and that 






g, e L\R^), 


£•0 > 0, go\B\n„ =0. 












Then 






0{T,-)\B\n, = 


- for any r € [0, T]. 













Proof: 

We use the level set approach to describe the interface Tr, see Osher and Fedkiw [16,. To this end, we introduce 
a function d = d{t, x) defined as the unique solution of the transport equation 

dtd + V^dit, x) • V = 0, i > 0, x e i?^ 

with the initial data 

r > for a; e B \ f7o, 
d{0,x) = do{x)(.C°°{R'), do{x)^\ _ , V^rfo 7^ on To. 

[ < for x e f^o U (i?3 \ B) 

Note that the interface F^ can be identified with a component of the level set {d(T, •) — 0}, while the sets B \flr 
correspond to {d{T, •) > 0}. Finally, 

\/^d{T,x) = X{T,x)n{x), for any x G F^, X{t,x) > for r e [0,r]. (4.13) 

For a given ^ > 0, we take 

as a (Lipschitz) test function in the variational formulation (|4.1ip to obtain 



min ^ —d;l 



/ gip{T,-) dx^ - / (gdtd+ gu-W^d) dx dt. (4.14) 

JB\nr S Jo J{0<d{t,x)<i} ^ ' 

Now, we have 

gdtd+ gu ■ V^d ~ g(dtd + u • V^dj = g{Y — u) ■ V^d 

where, by virtue of hypothesis (|4.12p and relation ()4.13|) . 

(V - u) • V^d e Wg '^(B \ fit) for a.a. t £ (0, r). (4.15) 

Introducing 

S{t, x) = distfl3 [x, d{B \ Vtt)] for t e [0, r], a; e B \ f^f , 

we deduce from (|4.15p and Hardy's inequality that 

i (V - u) • v^d e l2([o, t] X b \ g,). 





Finally, since V is regular, we have 

6{t,x) 



< c whenever < d{t, x) < ^; 



whence, letting ^ ^ in (|4.14p . we obtain the desired conclusion 



/ e(r,- 
JBxn^ 



- , . , dx = 0, 
10 



where we have used the fact that g G L°°{0, T; L'^{B)). 



Q.E.D. 



Thus, by virtue of Lemma [4. 1[ the momentum equation (|4.10p reduces to 

/ gu-ip{T,-) dx- (gu)o,5 • (^(0, ■) da; (4.16) 

^ f £iu-5fV5 + e[u(g)u] : Vj;V5+p(£i)div^<^ + (5gi'^div2,v3-/x f Vj;U + V^u- -diva-ulj : Vj;¥'j dx dt 

- / / A*cj ( V:eU + V^U - -diVa;Ul 1 : Vaj^J dx dt 

Jo JB\nt \ ^ J 

for any test function <p as in (j4.4p . We remark that it was exactly this step when we needed the extra pressure term 
6(f ensuring the density q to be square integrable. 

4.2 Vanishing viscosity limit 

In order to get rid of the last integral in (|4.16p . we take the viscosity coefficient 

/i = const > in Qt-, 
^^ -^ a.a. in ((0, T)x B)\ Qt- 
Denoting {guj,u^} the corresponding solution constructed in the previous section, we may use (|3.14p to deduce that 



Jilt 



while 



J B\nt 



VxU^ + V^u,^ - -diVj;Ut^I dx dt < c. 



V^^u^ + V^U(^ - -div^u^I dx dt < c, 



(4.17) 



where the latter estimates yields 



/ / fJ-u { ^xUuj + y^u^ - -diva;U;^I ) : V^(p dx dt = 
JQ JB\nt V -^ / 

/ / ^fJiZ^fiJZ, I Va;Ui^ + V^U(^ - -diva,u^I 1 : V^^tp dx di -^ as w ^ 
Jo JB\il.t \ •J / 

for any fixed y. 

On the other hand, as we know from Lemma 14.11 that the density q^ is supported by the "fluid" region Qr, we 
can still use p. lip . (|4.17p . together with Korn's inequality to obtain 



T 

J at 



Vj;U^P dx di < c. 



11 



Repeating step by step the arguments of the preceding section, we let w — > to obtain the momentum equation in 
the form 

/ gu ■ ip{T, ■) dx - (£»u)o,A-¥3(0, •) dec (4.18) 

= / [gu- dftp + g[vL ® u] : Vxf + _p(£i)div2.y> + Sg^div^'P — §(V2;u) : V x^ ) dx At 

Jo Jilt ^ ^ 

for any test function tp as in (|4.4p . Note that compactness of the density is now necessary only in the "fluid" part Qt 
so a possible loss of regularity of u^^ outside Qt is irrelevant. 

4.3 Vanishing artificial pressure 

The final step is standard, we let (5 — )■ in (|4.18p to get rid of the artificial pressure term 5g^ and to adjust the initial 
conditions, see [5l Chapter 6]. However, the momentum equation identity (|4.18p holds only for the class of functions 
specified in (|4.4p . The last step of the proof of Theorem 12.11 is therefore to show that the class of admissible test 
functions can be extended by density arguments. This will be shown in the following part. 

4.3.1 Extending the class of test functions 

Consider a test function ip € C^{[0,T] x R^:R^) such that 

ip{T, •) • n|r^ = for any t. (4-19) 

Our goal is to show the existence of an approximating sequence of functions (^„ belonging to the class specified in 
(14.41) and such that 



IIVnllwi.-((o,T)xS;fl3) < c, ¥>„ -^ (p, dt^Pn ^ 9t¥', and Vj,(p„ -^ V^V a.a. in Qt- (4.20) 

Combining (|4.20p with Lebesgue dominated convergence theorem we may infer that (p belongs to the class of admissible 
test functions for (|2.3p . 

In other words, we have to find a suitable solenoidal extension of the tangent vector field ipItt inside fir- Since 
Ft- is regular, there is an open neighborhood Ut of Ft- such that each point x S W,- admits a single closest point 
b^(x) e F^. We set 

h(T, x) = (^(t, bT-(x)) for all x e Ut- 

Finally, we define 

w(r,x) =h(r,x)+g(r,x), 

where 

g(T, x) = whenever x G F,-, 

and, taking the local coordinate system at x so that 63 coincides with x — bT-(x), we set 

g(r,x) = [0,0,g%T,x.)],dx,gHr,^) ^ -dx,h\T,^) - dx,h^{T,^)- 

We check that 

div2;w(r, •) = in Ur, w(r, ^Ir, = vi'^, Olr,- 

Furthermore, extending w(t, •) inside fi,-, we may use smoothness of (p and F^ to conclude that 

we W^'°°{Qt)- 
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As a matter of fact, a (smooth) extension of <p, solenoidal in the whole domain Qt, was constructed by Shifrin 
[I9l Theorem 4]. 
Writing 

<p = ((^ - w) + w, 

we check that the field w belongs to the class (|4.4p . while 

{(f - w)(t, •)|an^ = for any r > 0. 
Thus, finally, it is a matter of routine to construct a sequence a„ such that 

a„ e C,°°([0,r] X B; R^), supp[a„(T, •)] C fi. for any r e [0,T], 
in particular a„ belongs to the class (|4.4p . and 

l|a„||H'i.°o((o,T)xB;_R3) < c, a„ ^ (v3 - w), (9fa„ -> 5t(y - w), and V^an -^ V^itp - w) a.a. in Qt- 

Clearly, the sequence 

y„ = a„ + w 

complies with (|4.20p . 

We have completed the proof of Theorem 12.11 

5 Discussion 

The assumption on monotonicity of the pressure is not necessary, the same result can be obtained for a non-monotone 
pressure adopting the method developed in |4]. 

As already pointed out, the technicalities of Section [4. 3. II could be avoided by means of the construction of special 
test functions used in [7]. However, it would be interesting to show that the pressure is bounded in some L''{Qt), 
with q > 1, meaning that the estimate (I3.18P holds in Qt- 

As pointed out in the introduction, the general Navier slip conditions (jl.l[) are obtained introducing another 
boundary integral in the weak formulation, namely 



f f K{u-V)-ip dS:, dt 
"'0 JFt 



Taking k — k(x) as a singular parameter, we can deduce results for mixed type no-slip - (partial) slip boundary 
conditions prescribed on various components of Fj . 

Last but not least, the method can be extended to unbounded (exterior) domains with prescribed boundary 
conditions "at infinity". 
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